s 

333.955 
E30aibc 
2001 


MONTANA  STATE  LIBRARY 


3  0864  1001  9404  5 


AQUATIC  INVERTEBRATES  AND  HABITAT  OF  BILLMAN 

CREEK, 
PARK  COUNTY,  MONTANA 

July  2000 


A  BIOASSESSMENT  STUDY 


STATE  DOCUMFNTS  COLLECTION 
l''Y  <-  1  2003 

MONTANA  ?--.TE  LIBRA'^Y 

1515  '    ^:!i  AVE 
HFIFM-.    "•"    TANA  59620 


RECEIVED 

APR  0  4  20UZ 

DEQ/PPA 
Monitoring  &  Data  Manage,T,ent  Bure^ 


A  report  to 

The  Montana  Department  of  Environmental  Quality 
Helena,  Montana 


by 

Wease  Bollman 

Rhithron  Biological  Associates 

Missoula,  Montana 


December,  2001 


INTRODUCTION 

Aquatic  invertebrates  are  aptly  applied  to  bioassessment  since  they  are  known  to 
be  important  indicators  of  stream  ecosystem  health  (Hynes  1970).  Long  lives,  complex 
life  cycles  and  limited  mobility  mean  that  there  is  ample  time  for  the  benthic  community 
to  respond  to  cumulative  effects  of  environmental  perturbations. 

This  report  summarizes  data  collected  on  July  21,2000  from  2  sites  on  Billman 
Creek,  a  tributary  of  the  Yellowstone  River  in  Park  County,  Montana.  Aquatic 
invertebrate  assemblages  were  sampled  by  personnel  of  the  Montana  Department  of 
Environmental  Quality  (DEQ).  Study  sites  lie  within  the  Montana  Valley  and  Foothill 
Prairies  ecoregion  (Woods  et  al.  1999).  A  multimetric  approach  to  bioassessment  such  as 
the  one  applied  in  this  study  uses  attributes  of  the  assemblage  in  an  integrated  way  to 
measure  biotic  health.  A  stream  with  good  biotic  health  is  "  ...  a  balanced,  integrated, 
adaptive  system  having  the  full  range  of  elements  and  processes  that  are  expected  in  the 
region's  natural  environment. . ."  (Karr  and  Chu  1999).  The  approach  designed  by  Plafkin 
et  al.  (1989)  and  adapted  for  use  in  the  State  of  Montana  has  been  defined  as  "...  an  array 
of  measures  or  metrics  that  individually  provide  information  on  diverse  biological 
attributes,  and  when  integrated,  provide  an  overall  indication  of  biological  condition." 
(Barbour  et  al.  1995).  Community  attributes  that  can  contribute  meaningfully  to 
interpretation  of  benthic  data  include  assemblage  structure,  sensitivity  of  community 
members  to  stress  or  pollution,  and  functional  traits.  Each  metric  component  contributes 
an  independent  measure  of  the  biotic  integrity  of  a  stream  site;  combining  the 
components  into  a  total  score  reduces  variance  and  increases  precision  of  the  assessment 
(Fore  et  al.  1995).  Effectiveness  of  the  integrated  metrics  depends  on  the  applicability  of 
the  underlying  model,  which  rests  on  a  foundation  of  three  essential  elements  (BoUman 
1998).  The  first  of  these  is  an  appropriate  stratification  or  classification  of  stream  sites, 
typically,  by  ecoregion.  Second,  metrics  must  be  selected  based  upon  their  ability  to 
accurately  express  biological  condition.  Third,  an  adequate  assessment  of  habitat 
conditions  at  each  site  to  be  studied  is  needed  to  assist  in  the  interpretation  of  metric 
outcomes. 

Implicit  in  the  multimetric  method  and  its  associated  habitat  assessment  is  an 
assumption  of  correlative  relationships  between  habitat  parameters  and  the  biotic  metrics, 
in  the  absence  of  water  quality  impairment.  These  relationships  may  vary  regionally, 
requiring  an  examination  of  habitat  assessment  elements  and  biotic  metrics  and  a  test  of 
the  presumed  relationship  between  them.  Bollman  (1998)  has  recently  studied  the 
assemblages  of  the  Montana  Valleys  and  Foothill  Prairies  ecoregion,  and  has 
recommended  a  battery  of  metrics  applicable  to  the  montane  ecoregions  of  western 
Montana.  This  metric  battery  has  been  shown  to  be  sensitive  to  impairment,  related  to 
habitat  assessment  parameters,  and  consistent  over  replicated  samples. 

Habitat  assessment  enhances  the  interpretation  of  biological  data  (Barbour  and 
Stribling  1991),  because  there  is  generally  a  direct  response  of  the  biological  community 
to  habitat  degradation  in  the  absence  of  water  quality  impairment.  If  biotic  health  appears 
more  damaged  than  the  habitat  quality  would  predict,  water  pollution  by  metals,  other 
toxicants,  high  water  temperatures,  or  high  levels  of  organic  and/or  nutrient  pollution 
might  be  suspected.  On  the  other  hand,  an  "artificial"  elevation  of  biotic  condition  in  the 
presence  of  habitat  degradation  may  be  due  to  the  paradoxical  effect  of  mild  nutrient  or 
organic  enrichment  in  an  oligotrophic  setting. 


METHODS 

Aquatic  invertebrates  were  sampled  by  Montana  DEQ  personnel  on  July  2 1 ,  2000 
from  two  sites  on  Billman  Creek.  Sample  designations  and  site  locations  are  indicated  in 
Table  1 .  The  sampling  method  employed  was  that  recommended  in  the  Montana 
Department  of  Environmental  Quality  (DEQ)  Standard  Operating  Procedures  for  Aquatic 
Macroinvertebrate  Sampling  (Bukantis  1998).  In  addition  to  aquatic  invertebrate  sample 
collection,  habitat  quality  was  visually  evaluated  at  each  site  and  reported  by  means  of 
the  habitat  assessment  protocols  recommended  by  Bukantis  (1998)  for  streams  with 
glide/pool  prevalence. 

Evaluated  habitat  features  include  instream  conditions,  larger-scale  charmel 
conditions  including  flow  status,  streambank  condition,  and  extent  of  the  riparian  zone. 
Scores  were  assigned  in  the  field  to  each  habitat  measure,  and  these  scores  were  totaled 
and  compared  to  the  maximum  possible  score  to  give  an  overall  assessment  of  habitat. 

Aquatic  invertebrate  samples  and  associated  habitat  data  were  delivered  to 
Rhithron  Biological  Associates,  Missoula,  Montana,  for  laboratory  and  data  analyses. 
In  the  laboratory,  the  Montana  DEQ-recommended  sorting  method  was  used  to  obtain 
subsamples  of  at  least  300  organisms  from  each  sample,  when  possible.  Organisms  were 
identified  to  the  lowest  possible  taxonomic  levels  consistent  with  Montana  DEQ 
protocols. 

To  assess  aquatic  invertebrate  communities  in  this  study,  a  multimetric  index 
developed  in  previous  work  for  streams  of  western  Montana  ecoregions  (Bollman  1998) 
was  used.  Multimetric  indices  result  in  a  single  numeric  score,  which  integrates  the 
values  of  several  individual  indicators  of  biologic  health.  Each  metric  used  in  this  index 
was  tested  for  its  response  or  sensitivity  to  varying  degrees  of  human  influence. 
Correlations  have  been  demonstrated  between  the  metrics  and  various  symptoms  of 
human-caused  impairment  as  expressed  in  water  quality  parameters  or  instream, 
streambank  and  stream  reach  morphologic  features.  Metrics  were  screened  to  minimize 
variability  over  natural  environmental  gradients,  such  as  site  elevation  or  sampling 
season,  which  might  confound  interpretation  of  results  (Bollman  1998).  The  multimetric 
index  used  in  this  report  incorporates  multiple  attributes  of  the  sampled  assemblage  into 
an  integrated  score  that  accurately  describes  the  benthic  community  of  each  site  in  terms 
of  its  biologic  integrity.  In  addition  to  the  metrics  comprising  the  index,  other  metrics, 
which  have  been  shown  to  be  applicable  to  biomonitoring  in  other  regions  (Kleindl  1995, 
Patterson  1996,  Rossano  1995)  were  used  for  descriptive  interpretation  of  results.  These 
metrics  include  the  number  of  "dinger"  taxa,  long-lived  taxa  richness,  the  percent  of 
predatory  organisms,  and  others.  They  are  not  included  in  the  integrated  bioassessment 
score,  however,  since  their  performance  in  western  Montana  ecoregions  is  unknown. 
However,  the  relationship  of  these  metrics  to  habitat  conditions  is  intuitive  and 
reasonable. 

The  six  metrics  comprising  the  bioassessment  index  used  in  this  study  were 
selected  because  both  individually  and  as  an  integrated  metric  battery,  they  are  robust  at 
distinguishing  impaired  sites  from  relatively  unimpaired  sites  (Bollrhan  1998).  In 
addition,  they  are  relevant  to  the  kinds  of  impacts  that  are  present  in  the  Billman  Creek 
watershed.  They  have  been  demonstrated  to  be  more  variable  with  anthropogenic 


Table  1.  Sample  designations  and  sampling  sites.  Billman  Creek.  July  21,  2000. 


Sample 
designation 

GPS  location 
Latitude                 Longitude 

Station  1 

45°39'28"                  110°39'28" 

Station  2 

45°38'31"                 110°34'21" 

disturbance  than  with  natural  environmental  gradients  (BoUman  1998).  Each  of  the  six 
metrics  developed  and  tested  for  western  Montana  ecoregions  is  described  below. 


1.  Ephemeroptera  (mayfly)  taxa  richness.  The  number  of  mayfly  taxa 
declines  as  water  quality  diminishes.  Impairments  to  water  quality  which  have  been 
demonstrated  to  adversely  affect  the  ability  of  maj^ies  to  flourish  include  elevated 
water  temperatures,  heavy  metal  contamination,  increased  turbidity,  low  or  high  pH, 
elevated  specific  conductance  and  toxic  chemicals.  Few  mayfly  species  are  able  to 
tolerate  certain  disturbances  to  instream  habitat,  such  as  excessive  sediment 
deposition. 

2.  Plecoptera  (stonefly)  taxa  richness.  Stoneflies  are  particularly 
susceptible  to  impairments  that  affect  a  stream  on  a  reach-level  scale,  such  as  loss  of 
riparian  canopy,  streambank  instability,  channelization,  and  alteration  of 
morphological  features  such  as  pool  fi-equency  and  function,  riffle  development  and 
sinuosity.  Just  as  all  benthic  organisms,  they  are  also  susceptible  to  smaller  scale 
habitat  loss,  such  as  by  sediment  deposition,  loss  of  interstitial  spaces  between 
substrate  particles,  or  unstable  substrate. 

3.  Trichoptera  (caddisfly)  taxa  richness.  Caddisfly  taxa  richness  has  been 
shown  to  decline  when  sediment  deposition  affects  their  habitat.  In  addition,  the 
presence  of  certain  case-building  caddisflies  can  indicate  good  retention  of  woody  ? 
debris  and  lack  of  scouring  flow  conditions. 

4.  Number  of  sensitive  taxa.  Sensitive  taxa  are  generally  the  first  to 
disappear  as  anthropogenic  disturbances  increase.  The  list  of  sensitive  taxa  used  here 
includes  organisms  sensitive  to  a  wide  range  of  disturbances,  including  warmer  water 
temperatures,  organic  or  nutrient  pollution,  toxic  pollution,  sediment  deposition, 
substrate  instability  and  others.  Unimpaired  streams  of  western  Montana  typically 
support  at  least  four  sensitive  taxa  (Bollman  1998). 

5.  Percent  filter  feeders.  Filter-feeding  organisms  are  a  diverse  group;  they 
capture  small  particles  of  organic  matter,  or  organically  enriched  sediment  material, 
from  the  water  column  by  means  of  a  variety  of  adaptations,  such  as  silken  nets  or 
hairy  appendages.  In  forested  montane  streams,  filterers  are  expected  to  occur  in 
insignificant  numbers.  Their  abundance  increases  when  canopy  cover  is  lost  and 
when  water  temperatures  increase  and  the  accompanying  growth  of  filamentous  algae 
occurs.  Some  filtering  organisms,  specifically  the  Arctopsychid  caddisflies 
{Arctopsyche  spp.  and  Parapsyche  sp.)  build  silken  nets  v^th  large  mesh  sizes  that 
capture  small  organisms  such  as  chironomids  and  early-instar  mayflies.  Here  they  are 


considered  predators,  and,  in  this  study,  their  abundance  does  not  contribute  to  the 
percent  filter  feeders  metric. 

6.    Percent  tolerant  taxa.  Tolerant  taxa  are  ubiquitous  in  stream  sites,  but 
when  disturbance  increases,  their  abundance  increases  proportionately.  The  list  of 
taxa  used  here  includes  organisms  tolerant  of  a  wide  range  of  disturbances,  including 
warmer  water  temperatures,  organic  or  nutrient  pollution,  toxic  pollution,  sediment 
deposition,  substrate  instability  and  others. 

Scoring  criteria  for  each  of  the  six  metrics  are  presented  in  Table  2.  Metrics  differ 
in  their  possible  value  ranges  as  well  as  in  the  direction  the  values  move  as  biological 
conditions  change.  For  example,  Ephemeroptera  richness  values  may  range  fi-om  zero  to 
ten  taxa  or  higher.  Larger  values  generally  indicate  fevorable  biotic  conditions.  On  the 
other  hand,  the  percent  filterers  metric  may  range  fi-om  0%  to  1 00%;  in  this  case,  larger 
values  are  negative  indicators  of  biotic  health.  To  facilitate  scoring,  therefore,  metric 
values  were  transformed  into  a  single  scale.  The  range  of  each  metric  has  been  divided 
into  four  parts  and  assigned  a  point  score  between  zero  and  three.  A  score  of  three 
indicates  a  metric  value  similar  to  one  characteristic  of  a  non-impaired  condition.  A  score 
of  zero  indicates  strong  deviation  fi-om  non-impaired  condition  and  suggests  severe 
degradation  of  biotic  health.  Scores  for  each  metric  were  summed  to  give  an  overall 
score,  the  total  bioassessment  score,  for  each  site  in  each  sampling  event.  These  scores 
were  expressed  as  the  percent  of  the  maximum  possible  score,  which  is  18  for  this  metric 
battery. 

Table  2.  Metrics  and  scoring  criteria  for  bioassessment  of  streams  of  western  Montana 
ecoregions  (Bollman  1998). 


Score 

metric 

3 

2 

1 

0 

Ephemeroptera  taxa  richness 

>5 

5-4 

3-2 

<2 

Plecoptera  taxa  richness 

>3 

3-2 

1 

0 

Trichoptera  taxa  richness 

>4 

4-3 

2 

<2 

Sensitive  taxa  richness 

>3 

3-2 

1 

0 

Percent  filterers 

0-5 

5.01  -  10 

10.01-25 

>25 

Percent  tolerant  taxa 

0-5 

5.01  -10 

10.01-35 

>35 

The  total  bioassessment  score  for  each  site  was  expressed  in  terms  of  use-support. 
Criteria  for  use-support  designations  were  developed  by  Montana  DEQ  and  are  presented 
in  Table  3a.  Scores  were  also  translated  into  impairment  classifications  according  to 
criteria  outlined  in  Table  3  b. 

In  this  report,  certain  other  metrics  were  used  as  descriptors  of  the  benthic 
community  response  to  habitat  or  water  quality  but  were  not  incorporated  into  the 
bioassessment  metric  battery,  either  because  they  have  not  yet  been  tested  for  reliability 
in  streams  of  western  Montana,  or  because  results  of  such  testing  did  not  show  them  to  be 
robust  at  distinguishing  impairment,  or  because  they  did  not  meet  other  requirements  for 
inclusion  in  the  metric  battery.  These  metrics  and  their  use  in  predicting  the  causes  of 
impairment  or  in  describing  its  effects  on  the  biotic  community  are  described  below. 


The  modified  biotic  index.  This  metric  is  an  adaptation  of  the  HilsenhofFBiotic 
Index  (HBI,  Hilsenhoff  1987),  which  was  originally  designed  to  indicate  organic 
enrichment  of  waters.  Values  of  this  metric  are  lowest  in  least  impacted 
conditions.  Taxa  tolerant  to  saprobic  conditions  are  also  generally  tolerant  of 
warm  water,  fine  sediment  and  heavy  filamentous  algae  growth  (Bollman, 
unpublished  data).  Loss  of  canopy  cover  is  often  a  contributor  to  higher  biotic 
index  values.  The  taxa  values  used  in  this  report  are  modified  to  reflect  habitat 
and  water  quality  conditions  in  Montana  (Bukantis  1998).  Ordination  studies  of 
the  benthic  fauna  of  Montana's  foothill  prairie  streams  showed  that  there  is  a 
correlation  between  modified  biotic  index  values  and  water  temperature,  substrate 
embeddedness,  and  fine  sediment  (Bollman  1998).  In  a  study  of  reference 
streams,  the  average  value  of  the  modified  biotic  index  in  least-impaired  streams 
of  western  Montana  was  2.5  (Wisseman  1992). 

Taxa  richness.  This  metric  is  a  simple  count  of  the  number  of  unique  taxa  present 
in  a  sample.  Average  taxa  richness  in  samples  fi"om  reference  streams  in  western 
Montana  was  28  (Wisseman  1992).  Taxa  richness  is  an  expression  of  biodiversity, 
and  generally  decreases  with  degraded  habitat  or  diminished  water  quality. 
However,  taxa  richness  may  show  a  paradoxical  increase  when  mild  nutrient 
enrichment  occurs  in  previously  oligotrophic  waters,  so  this  metric  must  be 
interpreted  with  caution. 

Percent  predators.  Aquatic  invertebrate  predators  depend  on  a  reliable  source  of 
invertebrate  prey,  and  their  abundance  provides  a  measure  of  the  trophic 
complexity  supported  by  a  site.  Less  disturbed  sites  have  more  plentifiil  habitat 
niches  to  support  diverse  prey  species,  which  in  turn  support  abundant  predator 
species. 

Number  of  "dinger"  taxa.  So-called  "dinger"  taxa  have  physical  adaptations  that 
allow  them  to  cling  to  smooth  substrates  in  rapidly  flowing  water.  Aquatic 
invertebrate  "dingers"  are  sensitive  to  fine  sediments  that  fill  interstices  between 
substrate  particles  and  eliminate  habitat  complexity.  Animals  that  occupy  the 
hyporheic  zones  are  included  in  this  group  of  taxa.  Expected  "dinger"  taxa 
richness  in  unimpaired  streams  of  western  Montana  is  at  least  14  (Bollman, 
unpublished  data). 

Number  of  long-lived  taxa.  Long-lived  or  semivoltine  taxa  require  more  than  a 
year  to  completely  develop,  and  their  numbers  decline  when  habitat  and/or  water 
quality  conditions  are  unstable.  They  may  completely  disappear  if  channels  are 
dewatered  or  if  there  are  periodic  water  temperature  elevations  or  other 
interruptions  to  their  life  cycles.  Western  Montana  streams  with  stable  habitat 
conditions  are  expected  to  support  six  or  more  long-lived  taxa  (Bollman, 
unpublished  data). 


Table  3a.  Criteria  for  the  assignment  of  use-suppxjrt  classifications  /  standards  violation 
thresholds  (Bukantis,  1997). 

%  Comparability  to  reference 

Use  support 

>75 

25-75 
<25 

Full  support-standards  not  violated 

Partial  support-moderate  impairment- 
standards  violated 

Non-support-severe  impairment-standards 
violated 

Table  3b.  Criteria  for  the  assignment  of  impairment  classifications  (Plafkin  et  al.  1989). 

%  Comparability  to  reference 

Classification 

>83 

54-79 
21-50 
<17 

nonimpaired 
slightly  impaired 
moderately  impaired 
severely  impaired 

RESULTS 

Habitat  assessment 

Figure  1  compares  habitat  assessment  results  for  the  2  sites  visited.  Table  4 
itemizes  the  evaluated  habitat  parameters  and  shows  the  assigned  scores  for  each. 

Figure  1.  Total  habitat  assessment  scores  for  2  sites  on  Billman  Creek.  July  21,  2000. 


Station  1 


Station  2 


Table  4.  Stream  and  riparian  habitat  assessment.  Two  sites  on  Billman  Creek.  July  21, 
2000.  Incomplete  data  were  provided  for  Station  2;  the  score  for  that  site  is  based  only  on 
parameters  for  which  evaluations  were  available. 


Max.  possible 
score 

Parameter 

Station  1 

Station  2 

20 
20 
20 
20 
20 
20 
20 
20 
20 
20 

Bottom  substrate 
Pool  substrate  char. 
Pool  variability 
Channel  alteration 
Sediment  deposition 
Channel  sinuosity 
Channel  flow  status 
Bank  vegetation 
Bank  stability 
Vegetated  zone 

5 

6 

5 

20 

5 

20 

20 

10/10 

10/10 

8/8 

12 

11 

7 

9 

16 

6 

15 
8/8 
n.a. 
n.a. 

200 

Total 

137 

92 

Percent  of  maximum 
CONDITION* 

68.5 
SUB- 
OPTIMAL 

57.5 
SUB- 
OPTIMAL 

♦Condition  categories:  Optimal  >  80%  of  maximum  score;  Sub-optimal  75  -  56%;  Marginal  49  -  29%;  Poor 
<23%.   Adapted  from  Plafkin  et  al.  1988. 

Habitat  assessment  scores  indicate  sub-optimal  conditions  at  both  visited  sites. 
Instream  conditions  at  the  upstream  site  (Station  1),  however,  were  appraised  as  poor; 
monotonous  mud  substrate  was  reported,  with  pools  essentially  absent  due  to  fine 
sediment  deposition.  Streambank  stability  and  vegetative  protection  were  perceived  to  be 
optimal.  No  channel  alteration  was  reported. 

At  the  downstream  site  (Station  2),  instream  habitat  diversity  was  increased  by  the 
presence  of  filamentous  algae  and  root  mats,  though  much  of  the  substrate  was  reported 
to  be  mud.  Pool  variability  was  judged  marginal.  Extensive  channel  alteration,  in  the 
form  of  straightening,  was  reported,  thus  sinuosity  was  rated  marginal.  Information  about 
streambank  stability  and  riparian  zone  width  was  not  provided. 


Bioassessment 

Figure  2  summarizes  bioassessment  scores  for  aquatic  invertebrate  communities 
sampled  at  the  2  sites  in  this  study.  Table  5  itemizes  each  contributing  metric  and  shows 
individual  metric  scores  for  each  site.  Tables  3a  and  3b  show  criteria  for  impairment 
classifications  and  use-support  categories  recommended  by  Montana  DEQ. 

When  this  bioassessment  method  is  applied  to  these  data,  scores  indicate  that 
biotic  health  was  moderately  impaired  in  both  reaches  of  Billman  Creek.  The  upper  site 
did  not  support  designated  uses,  and  violated  standards,  while  the  lower  site  partially 
supported  standards. 


Figure  2.  Bioassessment  scores  for  2  sites  on  Billman  Creek.  July  21,  2000.  Revised 
bioassessment  metrics  and  criteria  (Bollman  1998)  used  as  reference. 


e 
E 
S 

M 
M 

s 


100 
90 
80 
70 
60 
50 
40 
30 
20 
10 
0 


Station  I 


Station  2 


Table  5.  Metric  values  and  bioassessments  for  2  sites  on  Billman  Creek.  July  21,  2000. 
Revised  bioassessment  metrics  and  criteria  (Bollman  1998)  used  as  reference. 


SITES 

Station  1 

Station  2 

METRICS 

METRIC  VALUES 

Ephemeroptera  richness 

1 

1 

Plecoptera  richness 

0 

0 

Trichoptera  richness 

2 

2 

Number  of  sensitive  taxa 

0 

0 

Percent  filterers 

7 

1 

Percent  tolerant  taxa 

19 

12 

METRIC  SCORES 

Ephemeroptera  richness 

0 

0 

Plecoptera  richness 

0 

0 

Trichoptera  richness 

1 

1 

Number  of  sensitive  taxa 

0 

0 

Percent  filterers 

2 

3 

Percent  tolerant  taxa 

1 

1 

TOTAL  SCORE  (max.=18) 

4 

5 

PERCENT  OF  MAX. 

22 

28 

Impairment  classification* 

MOD 

MOD 

USE  SUPPORT  t 

NON 

PART 

1 .  Classifications:  (NON)  non-impaired,  (SLI)  slightly  impaired,  (MOD)  moderately  impaired,  (SEV) 

severely  impaired.  See  Table  3b. 

'Use  suppoil  designations:  See  Table  3a. 


At  both  Billman  Creek  sites,  Ephemeroptera,  Plecoptera,  and  Trichoptera  taxa 
richness  were  all  much  lower  than  expected  for  a  foothill  stream.  No  sensitive  taxa  were 
present  in  the  sample.  The  proportion  of  tolerant  taxa  in  the  assemblage  was  greater  than 
expected. 

Aquatic  insect  communities 

The  modified  biotic  index  value  (6.27)  was  very  high  at  the  upstream  site  (Station 
1),  compared  to  expectations  for  streams  in  the  Montana  Valley  and  Foothill  Prairies 
ecoregion.  In  addition,  a  lone  individual  mayfly  was  present  in  the  sample,  and  no 
sensitive  taxa  were  collected.  These  findings  suggest  that  water  quality  was  severely 
impaired  at  this  site.  Impairment  could  be  due  to  nutrient  and/or  organic  enrichment, 
elevated  water  temperature,  or  both.  Habitat  deficiencies  appear  to  impact  benthic 
organisms  as  well,  the  caddisfly  fauna  was  limited  to  2  taxa,  and  seven  "dinger"  taxa 
were  collected.  Deposited  fine  sediments  apparently  compromised  instream  habitats.  The 
midge  Odonontomesa  sp.  was  relatively  abundant,  suggesting  that  a  portion  of  the 
benthic  substrate  was  sand.  No  stoneflies  were  present  in  the  sample.  This  could  indicate 
reach-scale  habitat  disturbances,  but  could  also  be  related  to  poor  water  quality,  for  which 
there  is  other  evidence.  The  presence  of  5  long-lived  taxa  suggests  that  periodic 
dewatering  or  other  catastrophes  have  not  impacted  the  site. 

At  the  downstream  site  (Station  2),  the  modified  biotic  index  value  (8. 15)  was 
even  higher  than  at  the  upper  site,  suggesting  that  water  quality  at  this  site  was  poor.  The 
single  mayfly  taxon  collected  here  was  the  tolerant  baetid  Centroptihim  sp.  Possible 
causes  of  the  impairment  were  the  same  as  those  suggested  for  the  upstream  site.  The 
dominant  organisms  here  were  immature  tubificid  worms,  suggesting  that  substrate 
elements  were  mainly  silt  and  mud,  and  lending  strength  to  the  hypothesis  that  water 
quality  severely  impaired  biotic  integrity.  Other  evidence  of  fine  sediment  deposition  was 
apparent  at  this  site  as  well;  only  2  caddisfly  taxa  and  6  "dinger"  taxa  were  collected. 
The  hemoglobin-bearing  midge  Cryptochironomus  sp.  was  common  at  the  site, 
suggesting  that  substrates  were  anoxic.  Three  long-lived  taxa  were  present  in  the  sample, 
and  none  were  abundant;  this  may  have  been  a  result  of  seasonal  dewatering. 

The  MT  DEQ  field  investigator  adds  the  following  comments:  "The  upper  site  is 
in  a  reach  that  has  a  silt  bottom  but  generally  appears  to  be  in  OK  to  good  shape. 
Nitrogen  loading  rates  appear  to  be  high  but  nuisance  algae  growths  weren't  present 
when  I  was  there.  The  lower  reach  is  degraded.  Sediments  are  anoxic,  nuisance 
filamentous  green  algae  growths  were  fairly  profijse  when  I  was  on  site,  and  it's  probably 
dry  periodically."  (Pat  Newby,  pers.comm.  Feb.  12,  2002). 

CONCLUSIONS 

•  The  performance  of  bioassessment  metrics  was  similar  at  both  sites  on  Billman 
Creek.  Severe  impairment  of  water  quality  is  suggested  by  the  results,  and  habitat 
inadequacies  were  apparent  as  well.  Water  quality  impairment  could  take  the  form 
of  nutrient  enrichment  or  temperature  elevation  or  both.  Habitat  disturbance 
included  fine  sediment  deposition. 

•  Periodic  dewatering  may  fiarther  compromise  biotic  integrity  at  the  downstream 
site  (Station  2). 


The  relationship  between  habitat  assessment  scores  and  bioassessment  scores  is 
illustrated  in  Figure  3.  The  red  curve  in  the  center  of  the  graph  represents  the 
hypothetical  relationship  between  habitat  quality  and  biotic  health  when  habitat 
degradation  is  the  sole  source  of  impairment  to  benthic  assemblage  health 
(Barbour  and  Stribling  1991).  Both  data  points  fall  below  the  line,  suggesting  that 
water  quality  was  a  major  component  of  the  impairment  evident  in  benthic 
assemblages  in  Billman  Creek. 


< 


Figure  3.  The  relationship  of  habitat  assessment  scores  and  bioassessment  scores  for  2 
sites  on  Billman  Creek,  July  21,  2000.  The  red  curve  represents  the  hypothetical 
relationship  between  habitat  scores  and  bioassessment  scores  if  habitat 
quality  solely  determined  biotic  health. 
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Taxonomic  data  and  summaries 

Billman  Creek 

July  21,  2000 
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Aquatic  Invertebrate  Taxonomic  Data 
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Site  Name:  Billman  Creek 
Site  ID:  Station  1  7/21/2000 
Taxon 


Approx.  percent  of  sample  lued:  20 
Quantity Percent 


HBI 


FFG 


Tubificidae  -  immature 

Sphaeriidae 

Physidae 

Acari 


Oreodytes  sp. 
Lara  avara 
Optioservus  sp, 

Zaitzevia  sp. 


10 

3.42 

9 

CG 

114 

39.04 

8 

CG 

34 

11.64 

8 

CG 

1 

0.34 

5 

PA 

Total  Misc.  Taxa 

159 

54.45 

Diphetor  hageni 

1 

0.34 

5 

CG 

Total  £pliemeroj)tera 

1 

0.34 

Brachycentrus  americanus 
Lepidostoma  sp.-turret  case  larvae 

9 
19 

3.08 
6.51 

1 
2 

OM 
SH 

Total  Trichoptera 

28 

9.59 

14 

4.79 

5 

PR 

4 

1.37 

4 

SH 

5 

1.71 

4 

SC 

1 

0.34 

4 

CG 

Total  Coleoptera 


24 


8.22 


Ceratopogoninae 
Tabanidae 
Hexatoma  sp. 


0.68 
0.34 

1.37 


PR 
PR 
PR 


Total  Diptera 


Cricotopus  Bicinctus  Or. 
Micropsectra  sp. 
Odontomesa  sp. 
Orthocladius  sp. 
Pagastia  sp. 
Parametriocnemus  sp. 
Polypedilum  sp. 
Rheocricotopus  sp. 
Radotanypus  sp. 
Tanytarsus  sp. 
Tvetenia  sp. 
Tribelos  sp. 


2.40 


Total  Chironomidae 


II 
3 

IS 
4 
3 
3 
3 
1 
2 

21 
1 
6 


3.77 
1.03 
5.14 
1.37 
1.03 

i:d3 

\M 
0.34 
0.68 

7.19 
0.34 
2.05 


7 

CG 

7 

CG 

4 

CG 

6 

CG 

1 

CG 

5 

CG 

6 

OM 

6 

OM 

4 

PR 

6 

CF 

5 

CG 

0 

PR 

73 


25.00 


Grand  Total 


292 


100.00 


Aquatic  Invertebrate  Summary  Data 


Site  Name:  Billman  Creel( 


Site  ID:  Station  1  7/21/2000 


TOTAL  ABUNDANCE 
Ephemeroptera  +  Plecoptera  + 
Trichoptera  (EPT)  abundance 

TOTAL  NUMBER  OF  TAXA 
Number  EPT  taxa 

TAXONOMIC  GROUP  COMPOSITION 


292 


26 
3 


GROUP 

Misc.  Taxa 

Odonata 

Ephemeroptera 

Plecoptera 

Hemiptera 

Megaloptera 

Trichoptera 

Lepidoptera 

Coleoptera 

Ehptera 

Chironomidae 


#TAXA       ABUNDAN  PERCENT 


4 

159 

54.45 

0 

0 

0.00 

1 

1 

0.34 

0 

0 

0.00 

0 

0 

0.00 

a 

0 

0.00 

2 

28 

9.59 

0 

0 

0.00 

4 

24 

8.22 

3,.. 

7 

2.40 

12 

73 

25.00 

CONTRIBUTION  OF  IX)MINANT  TAXA 

TAXON  ABUNDANCE 

Sphaeriidae 

Physidae 

Tanytarsus  sp. 

Lepidostoma  sp.-turret  case  lar 

Odontomesa  sp. 

SUBTOTAL  5  DOMINANTS 

Oreodytes  sp. 

Cricotopus  Bicinctus  Gr. 

Tubificidae  -  immature 

Brachycentrus  americanus 

Tiibelos  sp. 

TOTAL  DOMINANTS 


SAPROBIC  INDICES 
HilsenhofTBiotic  Index 


ICE  1 

PERCENT 

114 

39.04 

34 

11.64 

21 

7.19 

19 

6.51 

15 

5.14 

203 

69.52 

14 

4.79 

11 

3.77 

10 

3.42 

9 

3.08 

6 

2.05 

253 

86.64 

6.27 


RATIOS  OF  TAX  GROUP  ABUNDANCES 
EPT/Chironomidae 


0.40 


FUNCTIONAL  FEEDING  GROUP  (FFG)  COMPOSITION 
GROUP  #TAXA       ABUNDAN  PERCENT 


Predator 

6 

29 

9.93 

Parasite 

1 

1 

0.34 

Collector-gatherer 

12 

200 

68.49 

CoUector-filterer 

1 

21 

7.19 

Macrophyte-herbivore 

0 

0 

0.00 

Piercer-herbivore 

0 

0 

0.00 

Scraper 

1 

5 

1.71 

Shredder 

2 

23 

7.88 

Xylophage 

0 

0 

0.00 

Omnivore 

3 

13 

4.45 

Unknown 

0 

0 

0.00 

RATIOS  OF  FFG  ABUNDANCES 

Scraper/CoUector-filterer 

0.24 

Scraper/(Scraper  +  C.filterer) 

0.19 

Shredder/Total  organisms 

0.03 

DIVERSITY  MEASURES 
Shannon  H  (loge) 
Shannon  H  (log2) 
Evenness 
Simpson  D 


2.33 
3.36 
0.71 
0.18 


COMMUNITY  VOLTINISM  ANALYSIS 
TYPE  ABUNDANCE  PERCENT 

Multivoltine  57  19.35 

Univoltine  146  49.83 

Semivoltine  90  30.82 


#TAXA  ABUNDANCE  PERCENT 

Tolerant  5                     55           18.84 

Intolerant  1                       6            2.05 

Clinger  7                     54           18.49 
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Aquatic  Invertebrate  Taxonomic  Data 
Site  Name:  Billman  Creek 

Site  ID:  Station  2   7/21/2000  Approx.  percent  of  sample  used:  10 

Taion Quantity  Percent HBI         FFG 

Nematoda 

Tubificidae  -  immature 

Sphaeriidae 

Margaritifera  falcata 

Physidae 

Acari 


1 

0.32 

5 

PA 

232 

74.12 

9 

CG 

2 

0.64 

8 

CG 

1 

0.32 

4 

CF 

15 

4.79 

8 

CG 

7 

2.24 

5 

PA 

Total  Misc.  Taxa 

258 

82.43 

Centroptilum  sp. 

0.32 

2 

CG 

Total  Ephemeroptera 

0.32 

Corixidae  -  immature 

0.32 

8 

UN 

Total  Hemiptera 

0.32 

Brachycentrus  occidentalis 
Helicopsyche  borealis 

0.64 
0.32 

1 
7 

OM 
SC 

Total  Trichoptera 

0.96 

Dubiraphia  sp. 
Optioservus  sp. 

0.64 
0.32 

6 
4 

CG 
SC 

Total  Coleoptera 

0.96 

Tabanidae 

0.32 

8 

PR 

Total  Diptera 

0.32 

Cricotopus  Bicmctus  Gr.  2  0.64  7  CG 

Cryptochironomus  sp.  18  5.75  8  PR 

Heterotrissocladius  sp.  1  0.32  0  CG 

Klicrotendipes  sp.  1  0.32  6  CG 

Pagastia  sp.  13  4.15  1  CG 

Parametriocnemus  sp.  3  0.96  5  CG 

Polypedilum  sp.  5  1.60  6  OM 

Thienemannimyia  Gr. 3 0.96 6 PR 

Total  Cliironomidae 46 14.70 

Grand  Total  313  100.00 


Aquatic  Invertebrate  Summary  Data 
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Site  Name:  Billman  Creelc 


Site  ID:  Station  2  7/21/2000 


TOTAL  ABUNDANCE  313 

Ephemeroptera  +  Plecoptera  +                      ••  i 

Trichoptera  (EPT)  abundance                   ni  ^  4 

TOTAL  NUMBER  OF  TAXA                  v    J-  21 

Number  EPT  taxa  ',    3 

TAXONOMIC  GROUP  COMPOSITION 


GROUP 

#TAXA 

ABUNDAN  PERCENT 

Misc.  Taxa 

6 

258 

82.43 

Odonata 

0 

0 

0.00 

Ephemeroptera 

1 

1 

0.32 

Plecoptera 

0 

0 

0.00 

Hemiptera 

,  .».,.    '  - 

1 

1 

0.32 

Megaloptera 

-    '•-  t'TL^*^.-^  . 

Q 

0 

0.00 

Trichoptera 

0 

2 

3 

0.96 

Lepidoptera 

' 

0 

0 

0.00 

Coleoptera 

^. 

,  3 

0.96 

Diptera 

1 

1 

0.32 

Chironomidae 

8 

46 

14.70 

RATIOS  OF  TAX  GROUP  ABUNDANCES 
EPT/Chironomidae 

-  1    V  ,  «f 


0.09 


FUNCTIONAL  FEEDING  GROUP  (FFG)  COMPOSITION 


GROUP 


#TAXA       ABUNDAN  PERCENT 


Predator 

3 

22 

7.03 

Parasite                        — 

.„,:^ 

8 

2.56 

Collector-gatherer 

»™:ift>-*.-^ 

272 

86.90 

CoUector-filterer 

1 

1 

0.32 

Maci  ophyte-herbi  vore 

0 

0 

0.00 

Piercer-herbivore 

0 

0 

0.00 

Scraper 

2 

2 

0.64 

Shredder 

0 

0 

0.00 

Xylophage 

0 

0 

0.00 

Omnivore 

2 

7 

2.24 

Unknown 

1 

1 

0.32 

RATIOS  OF  FFG  ABUNDANCES 

Scraper/Collector-filterer 

2.00 

Scraper/(Scraper  +  C.filterer) 

0.67 

Shredder/Total  organisms 

0.00 

CONTRIBUTION  OF  DOMINANT  TAXA 

TAXON  ABUNDANCE 

Tubificidae  -  immature  232 

Cryptochironomus  sp. 

Physidae 

Pagastia  sp. 

Acari 

SUBTOTAL  5  DOMINANTS 

Polypedilum  sp. 

Parametriocnemus  sp. 

Thienemannimyia  Gr. 


TOTAL  DOMINANTS 


PERCENT 


18 

5.75 

15 

4.79 

13 

4.15 

7 

2.24 

285 

91.05 

5 

1.60 

3 

0.96 

3 

0.% 

2% 

94.57 

SAPROBIC  INDICES 
HilsenhofTBiotic  Index 


DIVERSITY  MEASURES 
Shannon  H  (loge) 
Shannon  H  (log2) 
Evenness 
Simpson  D 


8.15 


0.84 
1.21 
0.28 
0.44 


COMMUNITY  VOLTINISM  ANALYSIS 
TYPE  ABUNDANCE  PERCENT 

Multivoltine  43  13.82 

Univoltine  263  83.95 

Semivoltine  7  2.24 


IIITAXA  ABUNDANCE  PERCENT 

Tolerant  7                     39          12.46 

Intolerant  1                       1            0.32 

Clinger  6                     13            4.15 


